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We study the problem of heat transport in one-dimensional (ID) spin-chain systems weakly cou- 
pled to three-dimensional phonons and impurities. We consider the limit of fast spin excitations 
and slow phonons, applicable to a number of compounds of the cuprates family, such as Sr2CuC>3, 
where the superexchange J is much larger than the Debye energy, Qd- In this case the Umklapp 
scattering among the spin excitations is strongly suppressed for all relevant temperatures. We ar- 
gue that the leading scattering mechanism for the spin excitations at not too low temperatures is 
the "normal" (as opposed to the Umklapp) spin-phonon scattering in which the non-equilibrium 
momentum is transferred from the spin subsystem to phonons where it quickly relaxes through the 
"phonon bath" . Because of the lower dimensionality of the spin excitations it is only the momentum 
along the chains which is conserved in such a scattering. We find that this effect leads to a particu- 
lar momentum- and temperature-dependence of the spin-phonon relaxation rate valid for the broad 
class of low-dimensional spin systems. Subsequently we demonstrate that the spin-phonon relax- 
ation mechanism is insufficient for the low-energy, long-wavelength ID spin-chain excitations, which 
make the thermal conductivity diverge. We complete our consideration by taking into account the 
impurity scattering, which in ID cuts off the quasi-ballistic spin excitations and renders the thermal 
conductivity finite. Altogether, these effects yield the following spin-boson thermal conductivity 
behavior: k s oc T 2 at low temperatures, k s cx T _1 at intermediate temperatures, and k b — const 
at higher temperatures T ~ Ob. The saturation at higher temperatures is of rather non-trivial 
origin and we provide a detailed discussion of it. Our results compare very well with the existing 
experimental data for Sr2Cu03. Using our microscopic insight into the problem we propose further 
experiments and predict an unusual impurity concentration dependence for a number of quantities. 

PACS numbers: 75.10.Pq, 71.10.Pm, 72.10.Bg, 75.40.Gb 



I. INTRODUCTION 

Experiments. Recent experiments in low-dimensional 
quantum magneto 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 have 
revealed remarkably strong thermal transport anomalies 
associated with the low-dimensional spin degrees of free- 
dom. In particular, an anisotropic thermal conductiv- 
ity, comparable in magnitude to that of metallic sys- 
tems, was observed in quasi-one-dimensional (ID) spin- ~ 
chain and ladder compounds with a large part of the heat 
current attributed to magnetic excitations^ 2 * 3 ^* 5 * 6 * 7 - In 
the quasi-2D and some quasi- ID materials various other 
effects were observed including strong increase of the 
thermal conductivity by a modest magnetic field 8 910,11 
and enhanced scattering due to suppression of the gap 
in the spin-gaped materials** 2 ** 3 " The thermal trans- 
port anomalies in layered cuprates associated with 
the magnetic and/or stripe excitations have also been 

fOUndl 14 ! 15 ! 16 ! 17 ! 18 

Both the strength and the temperature range for 
these new effects are very different from their counter- 
parts in 3D magnetic materials, where the heat trans- 
port by well-defined magnetic excitations exist only 
well below the corresponding 3D ordering transition 
temperaturei**^*** 2 *^ Note that the thermal transport in 
low-dimensional magnets has been studied before*** 2 ** 2 *^ 
and the qualitative reasons for that difference has been 
understood. In 2D systems magnetic excitations can be 



very well defined in a paramagnetic phase because of the 
large in-plane magnetic correlation lengths and in ID 
systems the soliton-like excitations can exist without any 
long-range order m& 

Theories. However, only recently a considerable 
progress in the theoretical understanding of the trans- 
port properties of ID quantum spin systems has 
been made, 27 . 28 ! 29 ! 30 ! 31 . 32 ! 33 ! 34 ! 35 ! 36 ! 37 ! 38 The main fo- 
cus of these recent studies has been on the relation- 
ship of the spin transport and conservation laws; 27 ! 28 
specifically on the possibility of an ideal conducting 
state in one-dimensional integrable and non-integrable 
systems While this problem is of significant inter- 
est, the importance of the spin-phonon and spin-impurity 
couplings, which break down the integrability of the un- 
derlying spin-only models, has also been discussedi 3 ** 3 ** 3 *^ 
In particular, it was obtained within the memory ma- 
trix formalism that the interplay of Umklapp scatter- 
ing and spin-phonon coupling leads to the exponen- 
tial temperature dependence of the thermal conductivity 
k cx e T l T , with T* proportional to the phonon's Debye 
energy 9j>i 37 i 38 

Boltzmann equation approach. In the present work 3 *^ 
we consider the problem of anomalous heat transport 
in quasi- ID spin-chain systems under a somewhat dif- 
ferent prospective. First, we derive microscopically a 
model of ID spin-boson excitations interacting with 
the 3D phonon environment and impurities, where the 
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spin-boson representation of the ID Heisenberg model 
is obtained by performing the standard Jordan- Wigner 
transformation followed by bosonization. 40 It is assumed 
that the spin-boson velocity is large in comparison with 
the phonon velocity, v S> c. The limit v c 
(J 3> ®d) corresponds to the experimental situation 
in Sr2CuC>3, SrCuC>2, (La,Ca,Sr)i4Cu24C>4i, and other 
cuprate materials^i^ where J ~ 2000K and v/c ~ 10. 
This large difference in spin and phonon energy scales 
is also responsible for the fact that, experimentally, the 
"spin peak" in thermal conductivity occurs at the tem- 
peratures well above the "phonon peak" ^j^^^ With 
this assumption and in the limit of weak spin-lattice cou- 
pling we solve the Boltzmann equation for the spin-boson 
distribution function and find the spin-phonon relaxation 
time T sp . 

Normal vs Umklapp. Why r sp is important? Gener- 
ally, the relaxation of the heat current and the resulting 
finite thermal conductivity should be due to the Umk- 
lapp processes or any other processes which do not con- 
serve momentum. Since the characteristic "bandwidth" 
for the spin excitations in the real ID spin chains mate- 
rials of interest is very large (J ~ 2000K) in comparison 
with the experimental temperature range, the Umklapp 
scattering of spin excitations on themselves is strongly 
suppressed (as ~ e _J / T ). Instead, by considering the 
limit v/c S> 1, we propose that the leading relaxation 
mechanism at not too low temperatures is the two-stage, 
bottle-neck process of (i) transferring momentum from 
the spin system to phonons, and (ii) subsequent dissipa- 
tion of the phonon momentum via an Umklapp process 
or impurity scattering. The central idea is that the excess 
momentum of a spin-boson waits the longest time to get 
transfcred to phonons, but once it is transferred it relaxes 
quickly. In other words, the phonon relaxation time t pp 
due to phonon-phonon scattering is much shorter than 
the spin-phonon relaxation time r sp , i.e. t pp <C t sp . As a 
result, the relaxation rate of such a two-stage process is 
determined by the "normal" spin-phonon scattering rate 
T 2-st = i T sp + Tpp]' 1 ~ T~ p l . Tne inequality t pp < r sp 
can be justified qualitatively with the help of the fol- 
lowing argumentation. First, the spin-lattice coupling is 
weak. Second, since c -C v the phonon states have a sig- 
nificant thermal population at the temperatures where 
the spin contribution to the transport is substantial. 

Furthermore, there is also a kinematic argument in fa- 
vor of T pp -c T sp that results from our analysis. The 
leading contribution to the transport comes from the spin 
bosons with a small momentum k <C T/v. We find that 
the most effective scattering mechanism for such bosons 
is an absorption of a "thermal" phonon, that is a phonon 
whose energy is of the order of T and momentum is of the 
order Pt = T/c 3> T/v S> k. For finite T such phonons 
have finite relaxation time r pp . We will show that the 
spin bosons in question have a divergent relaxation time 
T sp (k) cx 1/k 2 . Therefore, for spin bosons with suffi- 
ciently small k's the relevant phonon relaxation time r pp 
is shorter than the spin-phonon relaxation time r sp not 



only because of the weak spin-phonon coupling, but also 
due to the smallness of the characteristic momentum of 
the spin-boson, Pt 3> k. 

We have performed a detailed study of the kinematics 
of the spin-phonon scattering processes and obtained the 
corresponding transport relaxation rate T~ p l ~ Ak 2 T 3 /v 3 
for T « 6j), and r" 1 ~ Ak 2 TQ 2 D /v 3 for T > & D , 
A is a constant related to the spin-phonon coupling, 
A cx (g sp /c) 2 . The characteristic power of k 2 in the re- 
laxation rate can be traced back to the fact that spin ex- 
citations are confined to lower dimensions than phonons. 
The latter implies that it is only the momentum along 
the chains which is conserved in the spin-phonon scat- 
tering. This, in turn, leads to a larger scattering space 
for phonons and makes relaxation rate of the spin bosons 
tx k 2 . Note that for the problem of the phonon-phonon 
scattering in 3D the phonon scattering rate is oc fc 4 be- 
cause the scattering space is restricted by the momentum 
conservation in all three dimensions^ The temperature 
dependence of the relaxation rate is due to the bosonic 
nature of both the spin excitations and the phonons. 
Thus, the spin-phonon relaxation rate obtained in this 
work should apply to other low-dimensional spin systems 
as well. 

Infrared divergence. However, one can find that the 
considered scattering mechanism becomes too weak at 
low energy and ineffective for the dissipation of the low- 
energy spin bosons. Namely, within the formalism of 
Boltzmann equation, oc k 2 leads to the infrared di- 
verging thermal conductivity k s oc Jdk/k 2 ~ k~ l \k , 
fco — > 0, the situation familiar from the phonon ther- 
mal conductivity in 3D insulators4ii 4 i 4 il& In that lat- 
ter 3D phonon problem the regularization of such a di- 
vergence is non-trivial and involves a consideration of 
the higher-order phonon processes^^i degeneracy of 
phonon branches at high symmetry points^ or scatter- 
ing on boundaries. 

Impurities. What physical effect can render k s finite 
in our case? In contrast with the 3D systems^ because 
of the ID nature of spins, even a weak impurity poten- 
tial will have a dramatic effect on the low-energy spin 
bosons. The impurities generate a relevant interaction 
in the renormalization group (RG) sense. That is, they 
scatter low-energy excitations very well destroying their 
quasi-ballistic propagation. Below the so-called Kane- 
Fisher scale Tkf the scattering is very strong and leads 
to the localization of ID excitations^ On the other hand, 
at the temperatures well above this scale T 3> Tkf the 
impurity scattering can be analyzed perturbatively. As 
a consequence of such an analysis we find that impuri- 
ties result in the momentum-independent scattering rate 
T imp ^ n-T^ 1 , n being the impurity concentration. Total 
scattering rate t/~ \ = r^ 1 + r~^ p is finite at k = 0, which 
makes k s finite. 

Results. Altogether, the combined effect of impuri- 
ties and the spin-phonon scattering leads to the follow- 
ing spin-boson thermal conductivity behavior: (i) in the 
low-temperature, impurity-scattering dominated regime 
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k s oc T 2 , (ii) in the spin-phonon scattering domi- 
nated regime, intermediate temperatures T rn <C T <C 
Qd, k-s oc T^ 1 , and (m,) in the spin-phonon scatter- 
ing dominated regime, higher temperatures T > @d, 
k s = const oc T m /0£). Here the temperature 
T m corresponds to the maximum in k s (T). It is also 
the crossover temperature between the impurity- and 
phonon-scattering dominated regimes. 

We find that this temperature behavior of the ther- 
mal conductivity agrees very well with the available ex- 
perimental data for the spin-chain material Sr2Cu03- 5,6 
The thermal conductivity for the zig-zag chain SrCuC>2 
and the spin-ladder compounds (Ca,La,Sr)i4Cu2404i 
whose underlying spin models are different from the ones 
considered in this work, is briefly discussed. 

We analyze the impurity concentration dependence of 
several quantities. As a result we predict an unusual 
behavior of: (i) the crossover temperature T m oc n 1 / 6 , 
(ii) the spin-boson thermal conductivity maximum value 
K max _ Ks (y m ) n~ 2 ' 3 , and (Hi) the asymptotic value 
of n s at T > T m : nf oc n" 1 / 2 . 

Outline. Our paper is organized as follows. In Sec. [H] 
we introduce the spin-chain Hamiltonian and derive the 
spin-phonon and spin-impurity interaction terms. In Sec. 
IHII we find the spin-boson relaxation time due to scatter- 
ing off the phonons. In Sec. II VI the impurity contribution 
to the spin-boson relaxation time is calculated. In Sec. 
Ivl the thermal conductivity vs temperature is obtained 
and the results are compared to experimental data. In 
Section E] we also put forward several theoretical predic- 
tions for the impurity concentration dependence of differ- 
ent quantities and suggest further experiments. In Sec. 
I VII we discuss the notion of the mean free path and its ap- 
plicability to the transport by the long-wavelength spin 
bosons. We conclude by Sec. IVIII which contains the dis- 
cussion of our approach and approximations. Technically 
involved details concerning the spin-phonon collision in- 
tegral and impurity scattering are described in Appendix 
lAland Appendix iBl respectively. 



II. SPIN-PHONON AND SPIN-IMPURITY 
INTERACTION HAMILTONIANS 

Spin- chain Hamiltonian. The Hamiltonian of a single 
Heisenberg antiferromagnetic spin-i chain is: 



ffchain — J ^ Si • Sj_|_i. 



(1) 



This Hamiltonian, after the Jordan- Wigner transforma- 
tion, can be expressed in terms of fermionic operators tpi 
as follows: 



It is convenient to bosonize this Hamiltonian. To do that 
one introduces chiral fermionic fields iPl,r(x) such that 
ip(x) w tpL{x)e ikFX + %l)R(x)e~ lkFX . These fields can be 
rewritten in terms of spin boson field and its dual 

field 6(2) as 



1 



L,R 



nil 



(3) 



Using $, one writes the spin-chain Hamiltonian as: 



H = - I dx 



(4) 



L /2 

where J dx = f_ L / 2 dx, L is the linear size of the system, 
the Luttinger-liquid parameter IC is equal to 1/2, and the 
spin-boson velocity v is given by 



v = — J a. 
2 



(5) 



where a is the chain lattice constant. Furthermore, the 
spin-boson variables are redefined according to the rules: 



e = /c 1/2 e, 
$ = /cr 1/2 $. 



(6) 



Then we introduce the creation and annihilation opera- 
tors for the field 



^/2L\k\ 



(b{ + b_ k ) 



which diagonalize the Hamiltonian Ho, Eq. (@J: 
Ho=vJ2\k\b\b k . 



(7) 



(8) 



For details of this procedure see Ref. |4(J We would like to 
note that the Hamiltonian for the spin-i chain contains 
other terms beyond the LL form above. These terms do 
not lead to the finite heat conductivity by themselves 
and thus are not included in our consideration. The role 
of such terms in the presence of phonons via a phonon 
phonon-assisted Umklapp process has been considered 
elsewhere i2£ 

Spin-Phonon Hamiltonian. Since the superexchange 
J is a function of inter-site separation, the lattice vibra- 
tions are able to modify it. This mechanism will, there- 
fore, couple phonons to the spin degrees of freedom. We 
can account for this coupling by expanding J in gradients 
of inter-site distance ad x u x . 



J(a + ad x u x ) = J(a) + 5Jd x u x + 



(9) 



Hjw = 

i ^ 

+~ (24^ - l) (2Vj +1 V i+1 



(2) 



- 1 



where u x is the displacement along the chain. This trans- 
lates into the following interaction Hamiltonian: 

ffint = 9 -f J dx (K (d^f+JC- 1 (d K $) 2 ) (d x u x ) ,(10) 
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where g sp = SJa is a spin-phonon coupling constant. The 
atomic displacement vector u can be expressed in terms 
of phonon creation and annihilation operators as follows: 



u( r )= * y e 



£pe (al 



re + a -P( 



(11) 



where wp< = q|P| and summation runs over the three- 
dimensional wave-vector P and over three polarizations 
of 3D phonons denoted by I. Here £ is the polarization 
vector of a phonon, rrii is the mass of the unit cell, N is 
the number of unit cells in the sample, R = (x, y, z) is 
the position vector. Sound velocities q are much smaller 
than the spin-boson velocity v: 



Q <C v 



(12) 



For a chain specified by y = 0, z = the interaction 
Hamiltonian, Eq. 1)100. can be written, using Eqs. I|ll|) 
and J7J), as: 

Hs » = - lJ 7§ £ Vc(P,k,k')(a Pi blb k ,+h.c.) , (13) 



kk'Pi 



with 



V e (P,k,k') 



Pnkk' 



\/8rriiLopekk' 



k'-k 



(14) 



where k and k! are the ID momenta of spin bosons and 
P is the 3D momentum of a phonon. Note that only 
component of the total momentum along the chain is 
conserved, which is explicitly given by the 5-symbol. The 
projection of the polarization vector on the x-axis (£pe) x 
is equal to: 



pi) 



o 



longitudinal, 

transverse, 

transverse, 



(15) 



where two answers for the transverse phonons correspond 
to two possible choices of polarization It is convenient 
to choose the first of these polarizations to lie in the plane 
given by x and P and the second one to be normal to this 
plane. Since the projection £ x is zero in the second case, 
such phonons do not couple to the spin bosons and are 
not discussed in this paper any further. 

Altogether, the Hamiltonian (|13fl describes emission 
and absorption of phonons by the spin-density excitations 
represented by the spin bosons, see Fig. ^ 

Spin- Impurity Hamiltonian. In addition to interact- 
ing with phonons spin bosons are also scattered on impu- 
rities. Microscopically, the origin of that effect is simple: 
impurity leads to a local variation of the superexchangc 
coupling J, which leads to scattering of magnetic excita- 
tions. The impurities affect the low temperature trans- 
port properties of the Tomonaga-Luttinger liquid in a 
dramatic way because they act as almost ideal backward 



scatterers for the low- lying excitations. The most impor- 
tant (relevant in RG sense) part of the impurity Hamil- 
tonian for an impurity located at xo can be written as: 

#imp = a6 J imp e~ tkFX0 ip{(x )i; R (x ) +h.c. (16) 



5J-, 



imp 



COS 



Although the effect of impurities is very strong for the 
low-lying excitations, for small SJi mp <C J and not too 
low temperatures the impurity scattering can be consid- 
ered perturbatively. 



III. PHONON MECHANISM OF RELAXATION 

Boltzmann equation. In order to study the relaxation 
of spin excitations on phonons we will solve the Boltz- 
mann equation for spin bosons coupled to the bath of 3D 
phonons. The stationary Boltzmann equation has the 
form: 



vd x f k = S k [f], 



(17) 



where f k is the spin-boson distribution function and 
Sk[f] is the collision integral. In general, the collision 
integral is a non-linear functional of both f k and the 
phonon distribution function npf. Sk — S k [f,n]. How- 
ever, as we argued in the Introduction, one can assume 
that the relaxation of phonons is quicker than that of 
the spin bosons. Thus, even in the presence of the 
temperature gradient d x T, phonons will be treated as 
if they are in a local equilibrium with themselves and 
this equilibrium is characterized by a local temperature 
T(x) = T + xd x T: 



npi 



np t (T(x)) 



1 



z Ui Pl /T{x) _ i ' 



(18) 



Therefore, one can write « Sk[f, n°) = £&[/]. 

Collision integral. The collision integral for spin 
bosons can be generally written as: 



S k = 



W w f k {f k , +1)- W k , k f k ,{f k + l) 



(19) 



where stands for J dk' /2ir and W kk ' is the total prob- 
ability of the spin excitation to be scattered from the 
state k to the state kl . For the processes of scattering 
due to phonons such probabilities are given by: 



Wkk- 



w 



kk 



,p(n°p e + l) 



w 



k'kP' L Pi 



(20) 



where the symbol J p stands for (l/2ir) 3 f d 3 P. Similar 
expression for the probability Wy k of scattering from 
kl to k can be obtained by permutation of k and k! 
in 1)20(1. Here the "elementary" scattering probabilities 
w kk'P °f t ne spin boson due to emission or absorption of 



the phonon with the momentum P and polarization £ are 
determined from Eqs. JT3JI, lfT^|l : 



w 



kk'P 



g 2 s P Vo 

8rrii 
x6(k' 



P^\kk' 



-ml 



P\\ - k)8{w k , 



(21) 



" Wp; — LVk ) 



where Vq is elementary cells volume, ujk = v\k\ is the 
spin-boson energy. In all the collisions the total energy 
and the projection of the momentum along the chain is 
preserved. These conservation laws are enforced by the 
(5-functions in Eq. l(2*T|) . 

Linearized Boltzmann equation. Boltzmann equation 
(|17|l should be solved in the presence of a non-zero tem- 
perature gradient d x T ^ 0. Assuming that the gradi- 
ent d x T is small one can introduce a usual ansatz for 
the distribution function fk to linearize the Boltzmann 
equation: 



fk — fk + fk ; 

fk(T) = 



(22) 



e ^ fc /r_l' 



where fl is a non-equilibrium correction to the equilib- 
rium distribution function Function fl is considered 
to be small. Since the collision integral vanishes identi- 
cally in an equilibrium state when fk = f and up = n p , 
one expects that fl should be proportional to the tem- 
perature gradient d x T. 

Thus, the Boltzmann equation for the spin bosons to 
the first order in d T T can be rewritten as: 



t { x > dk 



Sk, 



(23) 



where Sk now stands for the linearized collision integral 
which can be expressed using the notations of Eq. I|22|) 
as: 



S k = 



E 



s 



(1) 

kt 



s 



(1) _ 

kt ~ 



W hh , 



kk'P 



k'P 



s 



w 



k'kP 



k'P 



c(2; 



(n Pi + fl + 1) fl 
+ {fk - r&) fl 

+ f°k + 1) fl 

(fk' ~ T&t) fl 



(24) 



(25) 



(26) 



The integral S^ accounts for two types of collision 
events: (i) spin boson with the momentum k emits a 
phonon with the momentum P and scatters into the state 
with the momentum k'; (ii) process inverse to (i). Like- 
wise, S^ describes the absorption of a phonon by spin 
boson with the momentum k and the corresponding in- 
verse process, see Fig. ^ 





k' 



k' 



FIG. 1: The scattering diagrams of the spin boson on the 
phonon for the processes contributing to the collision integral 
(a) S^ 1 ' and (b) . Solid lines are spin bosons, wavy lines 
are phononsi— 



Once the Boltzmann equation (|23|l is solved for fl the 
spin-boson thermal current density JJ, can be written as: 



T s — 
j e — 



2 ; ji 

Z7T 



(27) 



where the coefficient k s is the spin-boson thermal con- 
ductivity. The total thermal conductivity observed in 
experiment is the sum of k s and the phonon thermal 
conductivity. 

Kinematic considerations. To solve the Boltzmann 
equation 112MI) one needs to evaluate integrals in Eqs. I|25() 
and H26(l . Part of this task can be accomplished without 
any approximations as the integral over P in S<-^ and 
S^ can be calculated explicitly. Because of the con- 
servation laws discussed above the remaining integration 
over k! will be restricted to some finite intervals whose 
extent and range depend on k. Mathematical details of 
these calculations are given in Appendix lAl 

The collision integral derivation can be simplified even 
more if we restrict our attention to the spin bosons whose 
energies are small: v\k\ <C min{T;0£>}. These energies 
are the most important since the transport is dominated 
by the bosons with small k. Such an approximation is 
used at various stages of Sk calculation. The consistency 
of our result with this latter assumption can be verified 
afterward. 

Relaxation time. One can write the collision integral 
in Eq. 1(2*4^1 as a sum of two terms: 



fk 



1 Sp 



(k) 



+ SS k 



(28) 



where the first term has the usual relaxation time form, 
while the second one does not. However, the second term 
can be estimated and it is shown to be small for small 
spin-boson energy 5Sk/Sk oc (vk/T) 2 (for details see Ap- 
pendix . 

Finally, with the help of these approximations, the 
collision integral can be written in the relaxation time 
form Sk ~ fl/ T sp(k), where T sp (k) is the transport relax- 
ation time. Since the transport is dominated by the spin 
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bosons with small momentum, one can use that small- 
ness to simplify significantly the corresponding expres- 
sion l|A14|l in the collision integral to obtain: 



1 



Avk 2 [ Qd/v dk' 



r sv (k) 



T 



(k 



l\3 



2tt 2sinli (vfc'/2T) 



(29) 



where A is a constant proportional to the spin-lattice 
coupling, sinh(...) comes from the bosonic distribution 
function, and powers of k and k' are from the scatter- 
ing amplitudes V(P, k, k') in Eq. (|14|l and due to the 
integration over P. 

What does this formula tell us about the most effective 
scattering process for a spin boson with small momen- 
tum? One can see that the major contribution to the 
integral in l|29[l comes from vk' ~ min{T; ©£>}. That is, 
the scattered spin boson is "thermalized" . Therefore, the 
most important scattering process is the absorption of a 
thermalized phonon. 

Further simplification of the problem can be achieved 
by exploiting the small parameter c/v <C 1. Because the 
spin bosons are fast, the energy and momentum conser- 
vation dictates that the majority of the phonons which 
interact with the spin subsystem must have their momen- 
tum almost normal to the chain direction. Indeed, the 
energy conservation gives that the phonon momentum 
P = ujp/c = (u/c)||fc| — where k and k! are the ID 
spin-boson momenta, while the momentum conservation 
along the chain gives |P|| | — \k — k'\. Therefore, the mo- 
mentum component perpendicular to the chain for a typ- 
ical phonon is: |Pj_| w P ~ (t>/c)|Pii| ^> \P\\\, except for 
the case of almost elastic backward scattering k! w —k. 
The latter scattering event does not contribute substan- 
tially to the transport relaxation time (see Appendix lA"|l . 
Since the momentum P is almost normal to the chain 
the corresponding projection of the polarization vector 
for the longitudinal phonons given in Eq. <|15|) is small: 
£ x i = 0(c/v). At the same time the transverse phonon 
polarization projection is £ x t = 0(1). Thus, the most 
effective spin-phonon scattering is due to the transverse 
phonons. 

We would like to rewrite Eq. i|29|) in the following form: 



] Ak 2 T 3 



(A-) 



v 



with the auxiliary function T(z) defined as: 

rZ x 3 dx 



r(z)=±p- h{z)=r 

Moo) J 
and A given by: 



2sinrr x/2 



(30) 



, (31) 



2tt m.j \ c / 



(32) 



Ji(cc) = 12£(3) where ((x) is the zeta-function. One can 
easily verify that: 



T(Q D /T) = 



1 



(e D /Tf 



for T <C ©d 
for T » Q D 



(33) 



where Qd — ®d/ y/ ~ 6d/4 plays the role of a 
crossover temperature. That is, for the temperatures well 
below the Debye energy T<C 6o the integration in Eq. 
(|30|l is restricted by sinh _2 (x/2) and the relaxation rate 
is given by: 



1 



AT 3 k 2 



T sp (k) 



(34) 



However, for the temperatures comparable to 9d and 
higher (in fact, for T > Sd/4) the integration limit in 
Eq. I|31l) becomes small. This leads to the change in 
the temperature dependence of the relaxation rate from 
cx T 3 to cx T: 



1 



AQ 2 D Tk 2 



bp(fc) 



(35) 



The relaxation rate for the spin boson due to scattering 
on phonons in the form given in Eq. I|30|l with the limit- 
ing cases given by Eqs. (|34|l and H35|l are the main results 
of this Section. 

It is instructive to inspect and compare our results with 
the well-known phonon-phonon scattering rate for the 3D 
phonons, studied in detail many years ago42*Sia£ For the 
3D phonon-phonon scattering rate one finds that it de- 
pends on the higher power of the wave- vector <x k , while 
its temperature dependence is the same. This is because 
the scattering space in phonon-phonon case is restricted 
by the momentum conservation in all three dimensions. 
In our case, it is only the momentum along the chains 
which is conserved in the spin-phonon scattering. This 
leads to fewer restrictions and a larger scattering space 
for phonons and makes the relaxation rate cx k 2 . Thus, 
the characteristic power of k 2 in the relaxation rate Eqs. 
(|29|I - H35|) can be traced back to the fact that spin exci- 
tations are confined to lower dimensions than phonons. 
This implies that our results for the spin-phonon relax- 
ation rate should be applicable to other low-dimensional 
spin systems as well. The temperature dependence of the 
relaxation rate is due to the bosonic nature of both the 
spin excitations and the phonons. 

Infrared divergence. With the collision integral in the 
relaxation time form the Boltzmann equation is trivially 
solved: 



x _ v\k\T sp (k) df k ^ T sp (fc) 



T 



Ok 



(d x T) , (36) 



where we used that (v\k\/T)df° /dk w — 1/k in the small 
k limit. 

Once fl is found the energy current density can be 
determined from Eq. I|27|l which yields the thermal con- 
ductivity: 

f 2 <^ dk v " [ T ' V dk (vt\ 

This expression diverges at small k giving rise to an 
infinite k s . The divergence happens because the scat- 
tering of spin bosons on phonons is not sufficiently 



strong for k — > to ensure the convergence of the low- 
energy contribution to the thermal conductivity. Such 
a situation is familiar from the phonon thermal con- 
ductivity in 3D insulators J n the 3D phonon 
problem the regularization of such a divergence is non- 
trivial and involves consideration of the higher-order 
phonon scattering processes^ 4 * 4 ^ degeneracy of different 
phonon branches at high symmetry points^ and inter- 
face boundary scattering. 

We would like to note here that such a divergence is 
not related to any conservation law or integrability of the 
problem. The scattering simply becomes too weak at low 
energies and is unable to equilibrate the excitations. A 
model example of such a behavior would be a system of 
free, noninteracting phonons in a continuum, which are 
scattered only by the point-like impurities according to 
the Rayleigh's law. 49 All the conservation laws are broken 
in this case. Yet, the heat conductivity is infinite in such 
a system since the Rayleigh scattering is ineffective at 
low energies, which leads to a diverging integral similar 
to our Eq. (|37jl . 

In our case, because of the ID nature of spin system 
the impurity scattering of spin bosons is very effective 
at low energies, as we will show in the next Section, so 
that even weak impurity potential renders the thermal 
conductivity finite. 



IV. IMPURITY SCATTERING 

In a ID system impurities have a dramatic effect on the 
low-energy excitations. The Hamiltonian for spin excita- 
tions in a ID spin chain interacting with a single defect 
is given in Eq. Since the impurity Hamiltonian 

is a relevant perturbation in the RG sense, defects scat- 
ters the low-energy excitations very effectively destroying 
their quasi-ballistic propagation. Therefore, the disorder 
in the magnetic coupling J will remove the infrared di- 
vergence in Eq. (|3"T|) . 

The impurity Hamiltonian in Eq. I|16|l is not a low- 
order polynomial of the bosonic field. Thus, it is im- 
possible to account for the impurity scattering within 
the formalism of the Boltzmann transport theory. In- 
stead, we will evaluate the spin boson life-time r; mp using 
the Green's function perturbative expansion in powers of 

There are two issues we must clarify before proceeding 
with this approach. First, below the Kane-Fisher temper- 
ature the impurity scattering is very strong and leads to 
the localization of ID excitations. 59 On the other hand, 
at the temperatures well above this scale the impurity 
contribution can be analyzed perturbatively. Thus, for 
the perturbation theory in powers of S J lmp to be valid 
the temperature must be bigger than the Kane-Fisher 
temperature^: 






T T 



FIG. 2: Second-order diagrams for the spin-boson scattering 
on impurities. Vertices represent interactions with impurities 
at times t' and r". Shaded ellipses emphasize that there is 
a summation over virtual states with any number of spin- 
boson excitations. This is due to the fact that the interaction 
Hamiltonian is exponential in the bosonic field. Formally, the 
ellipses correspond to ^i^H^') e i^i{o,r") ^ 



Second, it is important to note that, generally, the life- 
time of an excitation is not equivalent to the transport 
relaxation time. While the only processes which violate 
the conservation of the momentum contribute to the lat- 
ter, any kind of scattering shortens the former. However, 
since Eq. does not conserve momentum and since in 
any typical scattering event the spin-boson momentum 
changes drastically, there is no distinction between these 
two time scales in our case. This justifies the use of our 
approach. 

The correction to the single-boson Green's function 
from the impurity scattering is evaluated as follows. The 
lowest-order single impurity contribution to the Green's 
function is: 



fc (r)-2? ofc (r)« j ((bl(r) + b_ k (r)) 

r im P (T')H imp (T"){b kl (0)+bl k ,(0) 



(39) 



dr'dr" 



The right hand side of this equation is proportional to the 
self-energy. Since the impurity breaks the translational 
invariance the self-energy depends on two momenta. As 
usual, the translational invariance will be restored after 
averaging over the random impurity positions. Since the 
perturbation H- lmp , Eq. (|I6() . is an exponential in the 
bosonic field $ rather than a polynomial, the calcula- 
tion of the Matsubara average involves an effective sum- 
mation of an infinite number of terms, as schematically 
shown in Fig. Details of such a procedure are given in 
Appendix iBl As a result of our calculation the retarded 
self-energy is obtained though the analytical continuation 
of the Matsubara self-energy: 



n S J hnp^ 

2aJ\k\T 



(40) 



T » T KF = J 



imp 



"'-'imp 



(38) 



with n being the dimensionless impurity concentration. 
Note that this self-energy is pure imaginary. 

Once the self-energy is found one can obtain the 
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FIG. 3: The sketch of the spin-phonon and spin-impurity 
relaxation rates t^ 1 and r~ mp as a function of (a) temperature, 
(b) momentum. In (a) 6d = Sd/4. 



dressed Green's function 

D k<u> = 



D 



Ok,Ul 



1 - D Qk ^Y,x 



where 



D, 



2u k 



Ok,uj 



LOT - 0T 



(41) 



(42) 



From here it is possible to find the lifetime of the spin 
bosons by solving equations on the poles of the Green's 
function: 1 = D k^^ w - The yields the following ex- 
pression for Tj m p: 



' imp 



where 



A 2 oc nS J? 



imp ■ 



(43) 



(44) 



The details of the derivation are given in Appendix B. Wc 
would like to note that the impurity-induced relaxation 
rate tT cx T~ l has been first obtained by Oshikawa and 
Affleck, Ref. |52|, in the context of the theory of electron 
spin resonance in the S=l/2 quantum antifcrromagnetic 
chains. 

Unlike the spin-phonon scattering, the relaxation rate 
due to the disorder is independent of k. The impurity 
scattering provides an effective relaxation mechanism at 
low energies and thus removes the divergence of the ther- 
mal conductivity in Eq. (|3TJl by effectively cutting off the 
low-energy spin bosons. 

Our Fig. |31 gives a qualitative comparison of the tem- 
perature and k dependencies of the spin-phonon and 
spin-impurity relaxation rates r" 1 and t~[ . 



V. THERMAL CONDUCTIVITY 

Cumulative result. Finally, having at hands the trans- 
port relaxation rates of the spin-boson due to scattering 



on the phonons, Eq. I|30(l . and on impurities, Eq. I|43|) . 
we are adequately equipped to calculate the thermal con- 
ductivity of the spin chains. The total relaxation time is 
given by: 



1 



Ttot(fc) 



' imp 



' sp 



(k) 



(45) 



Then the solution of the Boltzmann equation which we 
had for the phonon scattering 1(36|1 should be modified 
to: 



,1 _ v\k\T tot (k) d/g 



As a result the thermal current is given by: 

- 3 f fclfcl dfl 



J s e= [ v 2 kfl = {d x T) V - f — 

Jk 1 Jk T irnr 



imp 1 b P 



(46) 



(47) 



Thus, using the explicit expressions for the spin-phonon 
and spin-impurity relaxation rates the following formula 
for the thermal conductivity can be obtained: 



vT 2 f J / T x 2 dx 



a(T) 



ttA 2 J 4sinh 2 (x/2) x 2 + a(T) 
where x = vk/T and 

„,5 A 2 

a(T) 



(48) 



(49) 



with r(z) defined in Eq. 1)31(1 . Since the temperature is 
always much smaller than J one can safely replace the 
upper limit in Eq. I|48l) by infinity. 

Our Eq. I|48() specifies a function with the following 
properties. Depending on the value of ©£> it either has 
a single maximum at some T = T m (larger 0_d) or is 
a monotonic function (smaller 0_d). In both cases this 
function vanishes as T 2 at T <C T m and saturates at 
T ^S> T m . These qualitative features of k s (T) can be eas- 
ily established. Without impurities (A = 0) the integral 
in Eq. (|48|l diverges at the lower limit. The impurity 
scattering cuts off this divergence. When the temper- 
ature is small the impurity scattering is much stronger 
than the typical spin-phonon scattering, which corre- 
sponds to a(T) ^> 1. Thus, n s oc T 2 at low temper- 
atures. In the limit of larger temperatures the impu- 
rity scattering is much weaker than the spin-phonon one 
and a(T) <C 1. If we assume that the temperature is 
still much less than the Debye energy T m < T « 6d 
then F(Od/T) = 1- Therefore, in this temperature 
range Eq. gHJ yields k s cx T 2 ^/a{T) cx 1/T. At yet 

higher temperatures T > ©d the temperature behavior 
of the spin-phonon relaxation rate changes from oc T 3 
to cx T, which is reflected in the change of r(©£>/T) 
from 1 to w (©1,/T) 2 , see Eqs. l|50 ) -ft5 |l . This gives 
k s cx T 2 y/a(T) cx T°. Altogether, assuming ©u is not 



too small: 



K a (T) oc 




for T«T m , 

for T m < T < & D , 

forT> e D . 



(50) 



where &d ~ ©d/4, as introduced in Sec. IIIII Note 
that the low-temperature, impurity-controlled thermal 
conductivity in the generic Luttinger liquids has been 
studied in Ref. [5^ using the memory-matrix formalism. 
Our oc T 2 result matches the result of that work for the 
Luttinger liquid parameter of the spin-i chain JC = 1/2. 
As one can see from Eq. H5()|l our results do not support 
suggestion of k s oc exp(T*/T) made in the experimental 
work, Ref. 0. Below we will discuss the origin of the 
maximum and of the saturation value of k s , but first we 
compare with the available experimental data. 

Comparison with experiments. To compare our results 
with experimental data it is convenient to rewrite the 
above expression for k s (T) , Eq. (|4"S1) , using the "reduced" 
temperature units t — T/T m with T m which corresponds 
to the maximum in the thermal conductivity: 



K s (t) 



7i 



dx 



1 



t 4 T(t)J 4sinh 2 (a;/2) x 2 +j 2 /t 6 T(t) 



,(51) 



where T(t) = T(^f-) from Eq. (j3T)) . and the constants 
71 and 72 are chosen in such a way that the temperature 
T = T m indeed corresponds to a maximum in k s and 
that K s (T m ) is indeed equal to n™ ax . For 6d not too 
close to T m the constants are: 71 = 1.26 and 72 = 1.83. 
Such a choice of variables simply replaces a combination 
of impurity- and lattice-related coupling constants en- 
coded in our A and A, whose actual values are generally 
unknown, by phenomcnological constants T m and K™ ax , 
found from experiment. In fact, in the limit ©£> — > 00 
this procedure completely determines our K s (t) given in 
Eq. (|51|l since there are no adjustable parameters left. 
The resulting K S (T) for this limiting case is shown in Fig. 
^]by the dashed line. 

Once we changed the variables to Eq. H51fl there is only 
one parameter to adjust: ®d- If we choose Qd — H.6T m 
we fit experimental data for the spin-chain compound 
Sr2Cu03 perfectly. For this material T m — 79AK and 
K max _ 36.7W m - 1 Fv~ 1 . Figure 0] shows the comparison 
of our results and experimental data, where the latter 
are shown by diamonds and the former is the solid line. 
Shaded area highlights the region where, experimentally, 
the phonon background subtraction creates a large un- 
certainty in the data£ 

Note that at high temperatures T > 6n Eq. i|51|) 
yields the asymptotic value of the spin-boson thermal 
conductivity 



7T7l 



e 



5.55 



D 



e. 



(52) 



Therefore, when we fix 0£> we fix nf as well. Our choice 



of Q D corresponds to k° 
marked by an arrow. 
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FIG. 4: Spin thermal conductivity k s normalized to its max- 
imal value vs reduced temperature T/T m . Diamonds are 
the experimental data for Sr2Cu03, Ref. 0- T m = 79. 4K 
and K™ ax = SQ.TWm' 1 ^ 1 for this material. Shaded area 
schematically represents the range where the phonon back- 
ground subtraction creates a large uncertainty in the data. 5 
Solid line is the results of this work, Eq. I51H . for Qd = 
11.6T m . Arrows mark the saturation value for the solid 
line at T > 9d and the spin-phonon scattering crossover 
scale 9u- Dashed line corresponds to n s (T), Eq. 151L for 
9d — > 00. [These results were recently presented in Ref. 1391] 



We would like to emphasize that apart from K™ ax and 
T m there is only one adjustable parameter, 65, which is 
used in our fitting procedure, and yet our theory yields 
an excellent agreement with the experimental data over 
the whole temperature range. 

Note that the transition between intermediate- and 
high-temperature regimes for k s is determined by the 
scale O d , which separates different temperature behavior 
for the spin-phonon scattering as discussed in Sec. IIIII 
It is this scale, not the Debye temperature itself, is a pa- 
rameter of our theory. For S^CuOa this scale is found 
to be Od « 3T m = 240K, shown in Fig. 01 by an arrow. 

The value of the Debye temperature &r> that comes 
out from our fit of thermal conductivity for Sr2Cu03 is 
Qd = H-6T m ~ 900K, which is higher than the esti- 
mations of Od for this material obtained by a different 
method. 5 Since in our theory the Debye temperature is 
a characteristic energy scale for the phonon bandwidth 
this may imply that the spin system is also coupled to 
the optical phonons. In fact, the typical width of the 
whole phonon spectra in the cuprates is about 900Ki£4 

To complete our comparison we provide here the 
asymptotic expressions for K S (T) in the low- and 
intermediate-temperature regimes. At small temperature 
Tkf C^C T m the thermal conductivity is: 



K S (T) ^ 2 7 i T 2 



372 T 2 n 



T 2 
2.26^3-. 

T 



(53) 



In this regime the conductivity is controlled entirely by 
impurities. 
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FIG. 5: Spin thermal conductivity k s normalized to its value 
at oo vs reduced temperature T/T* for various values of z = 
Qd/T* . (z decreases from the top to the bottom curve). 



When the temperature is higher the number of high- 
energy spin bosons increases. The most important relax- 
ation mechanism for these excitations is the spin-phonon 
scattering. However, the low-energy spin bosons are still 
scattered mostly by impurities. Thus, at intermediate 
temperature T m <C T <C Qd both mechanisms con- 
tribute to the relaxation of the thermal current: 



K S (T) 7T7i T„ 



2y^ T 



1.46- 



T 
~T~ 



(54) 



These asymptotic formulas can be tested experimen- 
tally, although the intermediate-temperature regime for 
Sr 2 Cu03 is not very well pronounced because 0_d is only 
3 times larger than T m . Such an intermediate asymptotic 
behavior might be more relevant for some other systems 
where the difference between T m and Qd is larger. 

For larger temperatures T > Qd the thermal conduc- 
tivity saturates. Thus, one of the most straightforward 
ways to verify our theory would be to measure the ther- 
mal conductivity for higher temperatures and check if the 
saturation really takes place as we predict. 

Different choices of Or). One may wonder at this point 
how much variation in the temperature dependence of k s 
can be expected from the theory given by Eq. (|48H or 
(|51|l if the parameters of the model, such as Qd, are var- 
ied. As is mentioned above there are three independent 
parameters: phonon- and impurity-couplings A and A, 
and the Debye temperature Qd- To generate a sequence 
of curves for the theoretical n s v.s. T the following pro- 
cedure is convenient. Since at T 3> Qd thermal conduc- 
tivity given by Eq. I|48fl saturates at a constant value, it 
is convenient to normalize k s (T) to its value at T = oo, 
i.e. nf, not to n™ ax . This is also because for Qd < 1.25 
k s (T) becomes a monotonic function and the determina- 
tion of n™ ax is ambiguous. Such a normalization reduces 
the number of free parameters to two. It is natural to 
choose one of those parameters to be a phonon crossover 



scale Qd (— ©d/4) and to define the other one as a refer- 
ence temperature T* . Using Eq. I|51|) with dimensionlcss 
constants 71 and 72 related through: 



1 = 



7T7l 



rp* 



2^71 Q d 



(55) 



one can generate a sequence of curves for n s /Kf v.s. 
t = T/T* for various values of z = Qd/T* shown in Fig. 
|S] Results for z =5, 4, 3, 2, 1.33, and 1 are shown (z 
decreases from the top to the bottom curve). The value 
of 72 = 1.83, same as in Fig. 0] is used to ensure that 
for Q D > T* (in reality, for Q D > 2T*) the reference 
temperature automatically corresponds to the tempera- 
ture of the maximum, i.e. T* = T m . Different choice 
of 72 would simply correspond to rescaling of T* leav- 
ing the form of K s (t) unchanged. Thus, Fig. [5] gives a 
description of possible variations in the results of the the- 
ory given by Eq. (|48|l within the experimentally relevant 
range of Qd/T* . 

As is described above, in order to fix the remaining 
two parameters and fit the experimental data one may 
associate the experimental value of T m with the position 
of the theoretical maximum and then vary Qd to match 
k™ for example. A similar strategy has been used 
to obtain Fig. 0] where the best fit closely corresponds to 
z = 3 curve in Fig. [3] If the experimental curve does not 
have a maximum^ one can choose any other reasonably 
distinct n s (T') in order to fix the choice of theoretical 
parameters T* and Qd- 

Origin ofT m . The temperature T m at which k s reaches 
its maximum can be expressed through the original pa- 
rameters of the problem as: 



Tn 



fv 5 A 2 \ 
V I2A) 



1/6 



1/6 



cx v 



vA 



(56) 



In a similar fashion one can relate k^' 
scopic parameters as: 



to the micro- 



vA 



2/3 



717T.4T4 



nS J 2 
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(57) 



The physical meaning of T m is the following. It is possible 
to define the "thermal" spin-phonon relaxation rate: 
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7r 



1 



T S p(k) 



AT 5 



(58) 



k — kq 



where fcr = T/v. By its definition tJL is independent 
of k. It characterizes rate of the spin-phonon collisions 
of the thermal spin bosons at a given temperature. The 
scale T m emerges naturally as a solution of the equation: 

T sp ~ r imp- (59) 

In other words, when T oc T m the thermal spin-phonon 
and spin-impurity collision rates are the same. At the 
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temperature higher (lower) that T m the spin-phonon (im- 
purity) scattering dominates. Note that in the low- 
temperature regime we assume that the temperature is 
still much larger than the Kane-Fisher temperature Tkf- 

Origin of . Our work predicts a remarkable be- 
havior of the thermal conductivity at the temperatures 
higher than the scale related to the Debye energy: satu- 
ration at a constant value, Eq. I|52|) . A parallel can be 
drawn with the thermal conductivity of metals, where the 
electronic part of k saturates at the temperatures above 
the Debye energy^ The parallel can be suspected to be 
even deeper because of the presence of the large energy 
scales in both problems, Fermi energy Ef ^S> T in the 
case of metals and J S> T in our case. 

A closer examination, however, reveals a difference be- 
tween the two cases. In metals, the thermal conduc- 
tivity temperature behavior is defined by the electronic 
specific heat C e (T) and the electron-phonon relaxation 
time T ep (T) through the quasiclassical relation: n e (T) = 
C e (T)vpT ep (T) /3. Note that it is the thermalized elec- 
trons which contribute most substantially to the trans- 
port. Specific heat is linear in T, while at T > Qjj the re- 
laxation time is inversely proportional to the temperature 
due to the thermal population of phonons, r ep (T) cx 1/T. 
Altogether, this renders n e = const. 

In our case the excitations are spin bosons, not 
fermions at large momentum hp. Therefore, in contrast, 
in the phonon-scattering dominated regime T ^> T m the 
spin-boson relaxation rate is strongly fc-dependent and 
the major contribution to the thermal current comes 
from the long- wavelength bosons with k <C T/v. Thus, 
the impurity-scattering continues to contribute to the 
temperature dependence of the thermal conductivity at 
higher temperatures. Therefore, in the case of spin chains 
the saturation of result of a non-trivial combina- 

tion of: particular temperature- and fc-dependencies of 
the impurity and phonon relaxation rates and the ID 
density of states of the spin-bosons. 

Further predictions. Our equations l|52|l . (|56l) and (|57|) 
allow us to formulate several predictions of our model for 
the impurity concentration dependence of several quan- 
tities that can be verified experimentally. 

First, we see that the temperature T m at which n s (T) 
reaches its maximum scales as n 1 / 6 . Such a weak depen- 
dence means that for a particular material the maximum 
in k s should be around the same temperature for a broad 
range of disorder concentration. 

Second, the maximum value of the thermal conductiv- 
ity K f ax = n s {T m ) scales as n" 2 / 3 , Eq. (JSZJl. Third, the 
saturation value has yet different concentration de- 
pendence n~ x / 2 . Such behaviors can be looked for in the 
materials with the isotope substitution for either mag- 
netic ions (Cu 2+ ) or surrounding ions as they all will 
lead to the local modification of the superexchange con- 
stant. If such dependencies are observed this will provide 
strong support to our theory. 

Another interesting experimental suggestion is based 
on the specifics of the spin-phonon scattering which we 



have found in Sec. II I II Namely, the most important scat- 
tering of spin bosons in the phonon-dominated regime is 
due to phonons whose momentum is almost normal to 
the direction of the spin chains. The characteristic en- 
ergy of such phonons is T. One can think of the following 
experimental setup: the heat current is directed along 
the chains while the thermal phonon pulse is induced 
perpendicular to the heat current (e.g., by a short laser 
pulse). This should lead to a sizable suppression of the 
spin-boson thermal current due to additional scattering 
by the extra phonons. The increase of the phonon heat 
current due to transverse pulse will be either negligible or 
can be accounted for by comparing to the results of the 
same experiment with the heat current and heat pulse 
both directed perpendicular to the chains. 

Consistency check. In considering the spin-boson re- 
laxation due to phonons we emphasized the importance 
of the two-stage, bottle-neck process, in which the spin 
boson scatters on the phonon via a "normal" process, 
while the momentum relaxation via an Umklapp process 
is taken care of by the phonon bath. Having carried out 
a comparison with experiments we can now verify the 
consistency of our results with the initial assumption. 
Namely, we can now check whether the characteristic 
spin-phonon relaxation time r sp is indeed much longer 
than the characteristic phonon-phonon relaxation time 

Tpp- 

Since below T m the scattering is impurity-dominated 
we need to estimate the relaxation times for T > T m 
only. We have already introduced the "thermal" relax- 
ation time for spin bosons rj p in Eq. (J5SJ. At T = T rn 
one can express such relaxation time through n™ ax , Eq. 
(|57jl . as: 



The phonon relaxation time can be estimated from the 
quasiclassical expression for the phonon thermal conduc- 
tivity: K p h = c 2 C(T m )T P p/3, where C(T m ) is the lattice 
specific heat and t pp is the average transport relaxation 
time of a phonon due to phonon-phonon scattering. 

Using these expressions we write the ratio of the ther- 
mal spin-phonon and the average phonon-phonon relax- 
ation times as: 




where the specific heat C(T m ) can be found as 
min{3; 127T 2 (T m / '9 D ) 3 / '5} i 3 . One can see that the small- 
ness of c/v in (|6*T|) is canceled by the largeness of &D/T m . 
For Sr2Cu03 in the region around T m the spin-boson 
and phonon parts of the thermal conductivity are of the 
same order. Also, C(T m ) 2. Taking the actual val- 
ues of «™ aa: , Kph, T m , J, etc. for S^CuOs we find that: 
T Tp/Tp P - 0.5-=- 1. 

It appears, superficially, that our assumption is not 
valid. However, at T > T m the typical momentum k 
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of a spin boson contributing to the thermal transport is 
much smaller than the "thermal" momentum kj- = T/v: 
k/kx ~ (T m /T) 3 . Thus, the typical relaxation time is 
much larger than the thermal one r sp ss t^- (kr/k) 2 . We 
also recall here that the phonons contributing to our two- 
stage process most effectively are the thermalized phonon 
with the momentum Pt ~ T/c as shown in Sec. II I II 
Therefore, r^ p of such a phonon is much shorter than the 
average transport relaxation time t pp which we estimated 
above, T pp -C t pp . Thus, the characteristic spin-phonon 
and phonon-phonon times involved in the two-stage pro- 
cess obey: 

» S « 1 , (62) 

due to both r sp 3> T^ p and T pp <C t pp . This proves the 
validity of our approach. 

We also estimated the Kane-Fisher temperature for 
Sr 2 Cu03 from 

nT KF 

which follows from Eqs. JSHJ, J^SJ, and (j57|) . Using 
experimental values for n™ ax and T m one obtains that 

^ « lO- 3 ^ 1 , (64) 

which means that our assumption for the impurity- 
controlled regime Tk f T m is valid unless the impurity 
concentration n is less than 0.1%. We are not aware of 
the actual level of disorder in S^CuOa, but such a high 
purity is unlikely. Therefore, this assumption is valid too. 

With the help of our asymptotic formula l|57|l one can 
estimate the spin-coupling constant using experimental 
values of K™ a:E and T m . This estimate leads to an un- 
reasonably high value of A/ J a 5 ~ 10 4 . First, a large 
value of the spin-phonon coupling constant can be due 
to an exponential dependence of the superexchange on 
the interatomic distance, as it is generated by the virtual 
tunnelling of electrons. Second, as one can see from Eqs. 
(|56|l and (|57|l the expression for the spin-phonon coupling 
contains high powers of both J and T rn . This creates an 
extra sensitivity of the result of our estimate to the un- 
known dimcnsionless constants, which our theory might 
be unable to captured 

Other materials. Other materials of the cuprate 
family in which the anomalous thermal transport 
has been observed include the zig-zag chain^ and 
spin-ladder materials, 1 and another spin-chain mate- 
rial, BaCu2Si20y, with much smaller superexchange 
constant X 

The most straightforward case for applying our theory 
seems to be the case of SrCu02 where the spin excitation 
spectrum remains gapless because of the frustrated cou- 
pling between the chains. Since the inter-chain coupling 
is of order J'/ J ~ 0.1 one may expect that while the low- 



and intermediate-T behavior of k s should be the same as 
in the case of Sr 2 Cu03, the high-T part might be differ- 
ent due to an additional scattering mechanism from the 
inter-chain interaction. However, as we show in Section 
IVII the high-T regime is fully dominated by the long- 
wavelength spin excitations, which means that J' should 
remain irrelevant and the behavior of k s (T) for SrCu02 
should be very similar to the "simple" chain material 
Sr2Cu03 at all temperatures. Nevertheless, we were un- 
able to fit the experimental data for SrCuO^ as success- 
fully as for Sr2CuC>3. On the other hand, as we learned 
from a very recent work Ref. 55 there might be an issue 
with the subtraction background for these materials. The 
thermal conductivity of the 5% Ca-doped SrCu02 does 
exhibit a saturation at high-temperatures, in an excellent 
agreement with our theory, but the peak at lower tem- 
peratures is much less pronounced. We therefore refrain 
from an explicit comparison of our theory to SrCu02 un- 
til the experimental issues are settled. 

In the case of (La,Ca,Sr)i4Cu2404i materials the ex- 
citations are gapped. Therefore, the results of our work 
cannot be straightforwardly applied to them because the 
gap, not the impurity scattering, controls the low-energy 
cut-off scale. Thus, one needs to reconsider the problem 
of the thermal transport for the gapped spin systems. On 
the other hand, the spin-phonon scattering part of such 
a consideration may stay very similar to our case, espe- 
cially for CagLa5Cu2404i. This latter compound pre- 
sumably contains a lot of 3D disorder (not related to the 
spin subsystem) due to a random distribution of Ca and 
La ionsji which might help to facilitate the fast phonon 
momentum relaxation and ensure that t pp <C t sp at high 
temperatures. We would like to note that another very 
recent experimental work, Ref. 1571 has discussed the role 
of the spin-phonon coupling in the gapped systems using 
a very different framework. 

The case of another spin-chain material, BaCu2Si207, 
is also different because of the lack of the spin and lat- 
tice energy scales separation. In this compound the spin 
boson velocity is of the same order as the sound veloc- 
ity, v ss cX Therefore, the kinematic considerations of 
our work will not be applicable for the collision integral 
in this system. In the temperature regime T ~ J, eas- 
ily achievable for this material, the thermal transport 
will be influenced by various other relaxation processes, 
which are negligible in our case. 

Altogether, the differences in the underlying spin mod- 
els for some of these materials will, most certainly, lead 
to different temperature behaviors of k s . Therefore, a de- 
tailed microscopic study in these systems along the lines 
of our work is necessary. 



VI. MEAN FREE PATH 

Having described successfully the temperature depen- 
dence of the thermal conductivity in Sr 2 Cu03 we would 
like to discuss the length scales that characterize scat- 



tering in this system. We also would like to address the 
question of whether the concept of the mean free path 
can be applied to describe thermal conductivity by spin 
bosons. More specifically, we would like to ask whether 
the thermal conductivity of spin bosons can be written 
in a simplified form (sometimes referred to as the simple 
kinetic equation): 



2500 



K s {T)=vC s (T)i(T) , 



(65) 



and if yes in what regime. Here C S (T) is the specific heat 
of the spin chains and £(T) is the mean free path. 

In the end, we also verify that the wavelength of a 
typical excitation contributing to the thermal current in 
Sr2CuC>3 is shorter than its mean free path. This con- 
dition is needed for the applicability of the Boltzmann 
equation formalism (quasiclassical approximation). 

The "total" relaxation time of the spin boson with the 
wave- vector fc is defined in Eq. I|45|) . Thus, one can in- 
troduce corresponding fc- and T-dependent length scale: 



£(k,T) = — T 



pen 



(k,T) 



(66) 



In order to be clear on terminology we would like to refer 
to this fc-dependent length scale as to the "fc-dependent 
mean free path" of the spin boson. At the same time the 
notion "mean free path" will be reserved for an averaged, 
only T-dependent, length scale. 

Cut-off length. Spin-phonon scattering rate in Eq. 166|) 
depends on k as t^ 1 oc k 2 while impurity scattering rate 
is fc-independent and plays the role of a cut-off scale. 
Therefore, it is natural to introduce the "cut-off" length 
as: 



£ imp (T) = £(0, T) 



vT 
A2 



(67) 



using explicit expression for rj^ , Eq. |@3J. It is inter- 
esting to find out the absolute value of such a length for 
the real system like Sr2Cu03. Since all the impurity- and 
phonon-related parameters, A, A, and ©£>, have been re- 
lated to the phenomenological constants T m (~ 80K) and 
K™ oa: (~ 37 W m _1 K _1 ), one can re-express A through 
them and obtain: 

7T 7l b C T T 

4n P (T) = — — j=- — ~ 300a — , (68) 



72 k 2 R T m V2 T n 



where a = 3.9A and b = 3.5A and c = 12. 7A are the 
lattice constants along and perpendicular to the chains 
in Sr2CuC>3, respectively^ Such a large value of £i mp 
(pa 1200 A at T = T m ) is in agreement with the esti- 
mates for the mean free path made in the experimental 
works A£ Note that the cut-off length grows as the tem- 
perature increases. This is a consequence of the impurity 
scattering in ID spin-chains which is non-perturbative in 
the number of spin bosons. 

k-dependent mean free path. We now express the fc- 
dependent mean free path (|66fl through the cut-off length 




FIG. 6: fc-dependent mean free path £(k,T) v.s. x — vk/T 
for three representative temperatures: T = T m (solid), T = 
T m /2 (dot-dashed), and T — 2T m (dashed). Arrows mark the 
cut-off length li mp (T) for each temperature. An approximate 
boundary of the bosonic thermal population is marked by the 
shaded area. 



(|67|l using the explicit form of the spin-phonon scattering 
rate, Eq. l|3Ujl: 



t(k,T) 



1 fvk\ 2 f T s " 



Tn 



72 V T 



D 



e 



,(69) 



where T(Qd/T) is the auxiliary crossover function intro- 
duced in Eq. I|31|) . which is related to the Debye tem- 
perature: r w 1 for T < D and T « (6d/T) 2 for 
T > 6j> Our Fig. shows the fc-dependent mean free 
path £(k,T) v.s. x = vk/T for three representative tem- 
peratures: T = T m (solid), T = T m /2 (dot-dashed), and 
T = 2T m (dashed). From Eq. ffity one can see that 
£{x, T) is a Lorentzian with the height given by the cut- 
off length ^imp(T) (marked by arrows in Fig. 0) and a 
temperature-dependent width. This width depends on 
T quite strongly: oc T~ 3 for T < 9_d and cx T~ 2 for 
T>Q D . 

The spin boson energy normalized by to the tempera- 
ture x = vk/T is a natural variable to describe thermal 
processes. The values of x which are of significance for 
such processes are < 3 because of the bosonic thermal 
population; this boundary is marked by the shaded area 
in Fig. |SJ One can see that the fc-dependent mean free 
path is essentially constant for x < 3 in the impurity- 
dominated regime T = T m /2, while it depends on x very 
strongly in the phonon-dominatcd regime T = 2T m . 

Impurity- dominated regime. Already from this obser- 
vation one can see that in the impurity-dominated regime 
all the spin bosons involved in the thermal transport are 
characterized by approximately the same mean free path 
£(k,T) « 4np(T), Fig. dash-dotted line, and thus the 
mean free path concept readily applies here. Namely, one 
can re- write our expression for the thermal conductivity, 
Eq. I]48p. somewhat differently, using the fc-dependent 
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FIG. 7: The temperature dependence of the cut-off length 
iimp (oc T, dashed) and the "thermal" length £th(T) — 
£(T/v,T) (dashed-dotted). In the phonon-dominated regime 
tth(T) cc T~ 5 for T < Q D and oc T~ 3 for T > Qd). The 
"effective" mean free path £ e //(T), Eq. 1751 . and the approx- 
imate "effective" mean free path £®^(T), Eq. 17611 . are also 
shown (solid and dotted, respectively). 



mean free path from Eqs. (|66[1 and (|69l) : 



T 



77 Jo (C* - 1> 



T) dx 



(70) 



Since £(a;, T) is ^-independent in the impurity-dominated 
regime, this expression can be reduced to the form 165|) 



k,(T) = v C S (T) £ imp (T), 



(71) 



with the mean free path equal to the cut-off length £(T) 
£im P {T), Eq. (JH8J, and the specific heat given by: 



C S {T) 



(e* - l) 2 



2 T 

3 J a 



,(72) 



where w = ttJ/2 and the well-known value of the integral 
(= 7r 2 /3) have been used. Thus, our result for the mean 
free path in the low-T regime (T < T m ) is £ oc T, which 
disagrees with experimental works where often a constant 
mean free path at low temperatures is reported^- 7 How- 
ever, this is the regime where experimental data for the 
spin part of the thermal conductivity are unreliable due 
to large phonon background. 

Phonon-dominated regime. One can see from Fig. [H| 
that the situation for £(x, T) is quite different in the 
phonon-dominated regime. While the cut-off length 
grows as the temperature increases, the range of the en- 
ergies of the spin bosons effectively participating in the 
thermal transport shrinks. In fact, it is a very unusual 
feature of the problem: as the temperature increases 
the thermal current becomes dominated by the excita- 
tions with longer and longer wavelengths. Therefore, our 
approximations related to spin excitations, such as the 
use of bosonization language, become better justified for 
higher temperatures. 



Using Eq. (|69|l one can define the "typical" momentum 
of the spin boson k, corresponding to the width of the 
Lorentzian £(k,T) — ^ mp /2: 



for T < 9, 



~ T T T 
Ja e D T 



J_ (In 

20a \ T 
60a T 



where the values J/T m — 20 and Gu/Tm = 3 for 
Si'2CuC>3 were used. 

Evidently, at T > T m one cannot possibly assume that 
all the spin-boson excitations are characterized by the 
same fc-dependent mean free path. Therefore, one can- 
not reduce Eq. (|70|l for k s (T) to a simplified form, Eq. 
(|65|l . In other words, the mean free path concept cannot 
be, strictly speaking, applied to the phonon-dominated 
regime of the problem. This is further demonstrated in 
our Fig. [7] which shows the temperature dependence of 
the cut-off length £- lmp (oc T, dashed) and the "thermal" 
length £th(T) = £(T/v 7 T) (fc-dependent mean free path 
for the excitation with the energy = T, dashed-dotted). 
The latter length, while indistinguishable from the cut-off 
length in the impurity-dominated regime T < T m , dies off 
very quickly in the phonon-dominated regime T > T m (as 
oc T- 5 for T < e D and then as oc T~ 3 for T > §d). This 
is simply to reiterate our point that the spin bosons are 
characterized by strongly fc-dependent mean free paths 
forT>T m . 

"Effective" mean free path. Nevertheless, one may 
want to insist on using the simplified equation (|65H in 
order to extract some "effective" mean free path £ e f f (T) 
characterizing the system, an approach used in the ex- 
perimental workSf*i*& 

Ka{T) = f v C S (T) £ e ff(T). (74) 

Using the explicit form of C s (T) obtained above H72|) the 
"effective" mean free path £ e ff(T) can be written as: 

3 



4//(T) 



clef 



(e x - If 



£(x, T) dx 



(75) 



This integral can be taken numerically and £ e ff(T) is 
plotted in Fig. [7] (solid line). Note again that all three 
lengths, £imp(T), £ t h(T), and £ e ff(T) are indistinguish- 
able for T < T m while they are very different for T > T m . 

Recalling that the height of the Lorentzian £(x, T) 
grows oc T while the width is decreasing oc T~ 3 for 
T < Qd and oc T~ 2 for T > Qd one can immediately 
conclude that £ e ff(T) should decay as oc T~ 2 for T < 6_d 
and oc T^ 1 for T > Qd as marked in Fig. Since 
C S (T) oc T, this yields n s oc const behavior for T > &rj- 

Using the fact that the width of £(x, T) is strongly In- 
dependent one can obtain an approximate expression for 
£ e ff(T) in Eq. I|75l) which is valid in both asymptotic 
limits T <C T m and T 3> T m : 



eff 



(T) = 



(T m ) 

Bt*y/T(f) 



rctan(Bt 3 v /f(i)), (76) 
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FIG. 8: The "effective" mean free path , Eq. (1751 . compared 
to the experimental results of Refs. I4l5l Asymptotic fits oc T 
for T < T m and oc l/T + C /T 2 for T > T m are shown by 
the dashed lines. Shaded area schematically represents the 
experimental error bars due to the phonon subtraction. 

where t = T/T m and constant B = 2.5. It is shown in 
Fig. □ (dotted line). Using explicit form of T(Q D /T) 
one can obtain all the correct asymptotic behaviors of 
4// (T) from it. 

Our Fig. [H] shows our "effective" mean free path, Eq. 
(|75|l . compared to the same quantity derived in Refs. Hlg 
from the experimentally measured k s (T) using Eq. (|74|l . 
Asymptotic fits oc T for T < T m and cx l/T + C/T 2 
for T > T m are also shown. Shaded area schematically 
represents the experimental error bars due to the phonon 
subtraction. 

Saturation in k s and parallel with metals. We would 
like to return to the discussion of the origin of saturation 
in n s (T) at T > Qd- One may be surprised, once again, 
by the seemingly identical features of the high-T thermal 
conductivity in both the spin chains and metals: (i) in 
both cases there is large scale (Ep and J, respectively) 
which yields "fast" excitations, (ii) in both cases specific 
heat oc T up to very high temperatures, (iii) for the satu- 
ration regime the dominant scattering is due to phonons 
in both cases, (iv) mean free path for metals and "effec- 
tive" mean free path for spin bosons are oc l/T in this 
regime. 

It is very tempting to conclude that in the case of the 
spin chains the l/T behavior in the "effective" mean free 
path is due to phonon population n p h oc T as in the case 
of metals. However, it is much more subtle in this case. 
As we discussed above and in previous Section, thermal 
transport in spin chains relies on the long-wavelength 
bosons, even more so for the high temperatures (see Fig. 
|SJ), and the T-dependence of the impurity scattering re- 
mains important. Simply put: if the impurity scattering 
would provide £ lmp v.s. T different from oc T, the In- 
dependence of the £ e ff would be different from oc l/T 
and there would be no saturation in k s (T). In fact, this 
seems to be the case for the related frustrated ladder 
(zig-zag chain) material SrCu02, which has gapless ex- 



citations as in spin chains, but its thermal conductivity 
does not seem to saturate at higher T£ 

Once again, it is an interesting and unusual situation 
in spin chains: it is always the long- wavelength range 
of excitations which is crucial for the thermal transport, 
temperature increase only enhances the phonon scatter- 
ing and shrinks the relevant energy region for spin bosons. 
Thus, our results are universal and do not depend on the 
manner in which bosonization or other approximations 
break down at higher energies. 

Quasiclassical approximation. It is instructive to verify 
the validity of the quasiclassical approximation which is 
implicitly used when utilizing Boltzmann equation for 
the lower-order scatterings. One needs to show for the 
"typical" spin boson that its wavelength A is shorter than 
its fc-dependent mean free path A < £. 

T > T m case. Consider T > T m first. We have 
already defined the "typical" momentum for this case in 
Eq. Using J/T m = 20 and G D /T m = 3 specific 

for Sr2Cu03 we find (neglecting numerical coefficients of 
order one): 

A~ a-L^^-^60a^- for T > Q D . (77) 

in m rn -*-m 

The "typical" ^-dependent mean free path is defined as: 
1= t(k, T) = £ J^I1 „ 150a-^. (78) 

Therefore, I > A is always fulfilled for the phonon- 
dominated regime. 

T < T m case. In the impurity-dominated regime 
"typical" spin boson has an energy v\k\ ~ T. Thus, the 
"typical" wavelength is: 

A~ a jL^i~20a^. (79) 
T m T T v ' 

The "typical" fc-dependent mean free path is simply a 
cut-off length: 

?=W T ) R> 300o^-. (80) 

Thus, £ = A at T& w T m /4 and the quasiclassical approx- 
imation is valid above this temperature. 

VII. CONCLUSIONS 

Approximations. We would like to review briefly the 
assumptions which have allowed us to solve the problem 
of spin-boson thermal transport in a reasonably simple 
analytical form. 

The main assumption is the smallness of the phonon 
relaxation time as compared to the spin-phonon scatter- 
ing time. This is justified since the spin-lattice coupling 
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is generically smaller than the phonon-phonon coupling. 
This also allows us to restrict ourselves with the low- 
est order in the spin-phonon coupling within the Boltz- 
mann equation formalism. In addition, in our case c <C v. 
Thus, at a given temperature the phonon population is 
much higher than that of the spin boson. This also im- 
plies higher relaxation rate for the phonons. If the above 
is true then the momentum relaxation of spin bosons 
can be viewed as a two-stage process with a bottle-neck. 
First, the momentum waits the longest time for a spin- 
boson— phonon collision to get transfered to the phonon 
subsystem. Second, the phonons quickly dissipate the 
transferred momentum. The first stage is the bottle-neck. 
Therefore, it controls the rate of the spin boson momen- 
tum dissipation. This is why we do not need to know 
microscopic details about the phonon-phonon collisions. 
It is enough to know that these collisions occur more of- 
ten than the spin-boson scattering on the phonons. 

The factor which might lead to violation of our as- 
sumption is exponential weakness of the phonon Umk- 
lapp relaxation at low enough temperatures. The Umk- 
lapp relaxation rate contains a characteristic exponent 
e -b&o/T ^ j g a con stant f or der of unity4£ Therefore, 
our approximation breaks down at very low temperature. 
Roughly speaking, this temperature corresponds to the 
point where the smallness of the spin-lattice coupling is 
comparable to the smallness of e~ be °/ T . In such a situ- 
ation the bottle-neck scheme described above will cease 
to work and it will be necessary to include a detailed 
description of the phonon-phonon collisions as well. In 
practice, however, the impurity scattering will become 
relevant at much higher temperatures. Our theory shows 
a very weak impurity concentration dependence of the 
temperature below which impurities control the trans- 
port: T m cx n 1 / 6 . This implies that even a very weak 
disorder will result in a substantial T m . 

We have verified our assumptions using the phe- 
nomenological parameters from the spin-chain material 
Sr2CuC>3 and obtained that the two-stage process ap- 
proximation is very well justified, see Sec. Ivl 

We would like to note that the problem of the spin 
transport in spin chains coupled to 3D phonons has been 
studied within a different approaches and gives a dif- 
ferent temperature dependence for the spin thermal con- 
ductivity: k s cx exp(T*/T). The essential difference be- 
tween our approaches is our treatment of the phonon- 
bath as a source of the momentum dissipation. If this 
process is neglected the "Normal" transfer of the momen- 
tum from spins to phonons does not lead to degrading of 
the thermal current. Then, it is only the "direct" spin- 
phonon Umklapp processes which will be contributing to 
the transport relaxation rate with a characteristic expo- 
nential dependence on temperature. However, this could 
only work if the phonon system would be very reluctant 
to dissipate the momentum, which can only happen at 
the temperatures of the order of the "phonon peak" ( 30K 
for Sr2Cu03), well below the temperature of the maxi- 



mum in the spin part of k (T m = 80K). 

Final remarks. In this paper we have studied the 
problem of anomalous heat transport in quasi-lD spin- 
chain systems coupled to the 3D phonon environment in 
the presence of weak disorder. We have derived a mi- 
croscopic model of ID bosonic spin excitations interact- 
ing with phonons and impurities using bosonization ap- 
proach. We have considered the spin-phonon scattering 
within the Boltzmann equation for the spin-boson distri- 
bution function. This equation has been solved in the 
limit of weak spin-lattice coupling and fast spin-boson 
excitations. Assuming that the phonon relaxation time 
T pp due to phonon-phonon scattering is much shorter 
than the spin-phonon relaxation time r sp we have ob- 
tained the transport relaxation rate for the spin bosons 

T tr l ~ ^p 1 <* k2 T 3 for T < and « k2T for T>B D . 
We have shown that within our model the spin-phonon 
scattering alone is insufficient to render the thermal con- 
ductivity finite because the spin-phonon relaxation mech- 
anism becomes ineffective for the spin bosons at low en- 
ergies. On the other hand, the ID impurity scattering 
provides a natural cut-off scale for the low-energy spin 
bosons and results in the momentum-independent scat- 
tering rate r^ p cx T _1 . Thus, the impurity scattering 
removes the divergence of the conductivity. We have cal- 
culated the thermal conductivity as a function of temper- 
ature and have shown that the low temperature transport 
is dominated by the impurity scattering while the high 
temperature transport is limited by both the impurity 
scattering and the spin-phonon collisions. Our main re- 
sults are in a very good quantitative agreement with the 
available experimental data. 

In implementing our approach we have also obtained 
an insight into various microscopic details of the problem. 
This has allowed us to formulate several predictions and 
suggest future experiments as well, see Sec. [V] One of 
the predictions is the saturation of the spin-boson ther- 
mal conductivity at high temperatures: a non-trivial re- 
sult due to impurity and spin-phonon scatterings and ID 
nature of spin-bosons. 

Further studies in the thermal conductivity of other 
spin systems using our approach are anticipated. 
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APPENDIX A: SPIN-PHONON COLLISION INTEGRAL 



In this Appendix we will present a detailed derivation of the relaxation time approximation for the spin boson 
collision integral. 

The integral over P in S^, Eq. can be evaluated explicitly. Let us denote the integrands of i|25|) other than 

delta-functions by h(P\\,u)p). Then, the integration over P can be performed giving the following answer: 

d(uj k , +UJ-P- u k )8{k' + P|, - k)h(P {] , up) = /i(fc - k', w k - u> k ,)6{\k\ - \k'\)6(x )e(e 2 D - so), (Al) 

with x^=v 2 {\k\-\k'\f -c 2 {k^k') 2 , (A2) 

where 0(x) is a step function and Qd is the Debye temperature. 

It is possible to act in the similar fashion for as well. The result of P integration for this part of the collision 
integral differs by permutation of k and k' . 

The step-functions 6 in IjAljl impose restrictions on the possible scattering states of the spin boson with the given 
momentum k. These restrictions are consequences of the conservation laws. Because of this step-functions the 
integration variable k' in Eqs. (|25|l and for the collision integrals SL' is bound to an interval C\ ' smaller than 
the Brillouin zone. This interval, in general, depends on k. To find ci we must solve: 

< x < Q 2 m , (A3) 

for |fc| > \k'\. Graphical solution of these inequalities for k > is presented in Fig. [5] as a lightly shaded area on 
(fc, k') plane. In this Figure lines OA and OB are bisectors k = ±fc'. Although not shown in Fig. |5J k < portion of 
the graphical solution can be obtained by an inversion of Fig. |5]with respect to the point O. 

Equivalently, we can express the solution as the following choice of the integration interval C» for \k\ > Qd/v: 



c (i) _j[k,k- Q Di /v] U [-fc + y/4( Q /i;) 2 fc 2 + (e De /v) 2 , -(1 - 2c t /v)k] for k > 0, . . 

1 \ [-(1 - 2c e /v)k, -k - ^A(c t /v) 2 k 2 + (Q m / V )2] [j[k + Q m /v, k] for k < 0, 1 ' 

where we use [k m i n , k max ] to denote the interval from fc m m to k max . For \k\ < Qdi/v: 

(i) _ / [fc, -(1 - 2c e /v)k] for k > 0, 



° e _ ( [-'(1 - 2c e /v)k, k] for k < 0. ^ 

(2) (2) 

For the contour C\ in the integral for S ki . 

c (2) \ [k + Q m/v, k] U [-(1 + 2c e /v) k, -k - @ m /v] for k > 0, 

1 \ [k - ^4(c e /v) 2 k 2 + (e De /v) 2 , (1 + 2c e /v)k] U[k,k- Qdi/v] for k < 0. ( } 

The fc > part of this solution corresponds to the darker shaded area in Fig. 

Keeping in mind the above results we rewrite the collision integral in the following form: 

s k = J2 s w> ( A7 ) 
kk'(k-ky 



s - 4 / dk ' 

C, ^ 



e (|fe|-|fe'|)/fe T _ i 



2 j e \k\/k T _ I i e -\k'\/k T - , | . 




where 



* = ^Lr (¥) , (A9) 



647r 3 mi \ q 

and k T = T/v . (A10) 



The contour Ci is the union U C^ 2; . For |fc| < Ooe/v we must integrate between ±(|fc| + ®Dt/v). For |fc| > Ooe/v 
we must integrate between k ± Qdi/v and between — fc ± ^/4(Q/u) 2 fc 2 + {Qdi/v) 2 - In both cases a window of the 
width 4cf|fc|/i> centered around fc' = — fc must be deleted from the integration interval. In Fig. |5|one can identify this 
window with the unshaded area around line OB. 
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B 



FIG. 9: This diagram shows the regions in the momentum space, allowed by the momentum and energy conservation, for the 
spin boson scattered from the state k to the state k' . The processes in which the original spin boson decays into another spin 
boson and emits the phonon are represented by the light gray area. For them \k'\ < \k\ as dictated by the energy conservation. 
The dark gray area corresponds to the processes of phonon absorption by the spin boson. For them \k'\ > \k\. Lines OA and 
OB are bisectors k' = ±fc. 



We start by making the following observation. In general, the intervals Cg in (| A8|) differ for different polarizations 
I = l,t: the size of the gap around k' = —k depends on the value of eg, so do @di/v which specify the outer limits 
of the integration intervals. However, we will show later that these variations in Cg do not change the integrals (|A8|) 
significantly. Therefore, it is convenient to introduce Cq = Ci n Ct and SCg = Cg — Co- Using this notation we write: 



S,= 



sc, 



(All) 



As it was mentioned the integrals over 5Cg in the above expressions are small. Hence, it is permissible to integrate 
over Co in both St and Si. In such a situation the summation over polarization I can be done easily since in both St 
and Si only the product Ag(£g) 2 depends on polarization. Thus: 



A = J2 A * (&)l = A u 



(A12) 



where £i x = c{\k — k'\/v\\k\ — |fc'|| and £ tx = a/I — (c 4 |fc — £/|/v||/c| — |fc'||) 2 . Therefore, the collision integral itself is 
equal to: 



Co 27r 



kk'{k-k') 2 



,(\k\-\k'\)/k T _ i 



,\k\/k T _ 1 i e -\k'\/k T _ 1 i ' 



z \k>\/k T _ r 



The first term of this integral can be calculated for small \k\: 



All I ^ 



Co 



2n 



kk'{k - k 



>Y2 



,(\k\-\k'\)/k T _ I 



,\k\/k T _ 1 l2 f e D /v dk , 

Ajk ;„ 



,\k'\/ki 



1 



'VI 3 



kr J~e D /v 27r 2cosh(fc'//c T ) 



which gives the spin-boson transport relaxation time: 



v 3 /AT 3 k 2 for T < B D 
~ | v 3 /AS 2 D Tk 2 for T > e D 



with A = hMj « 2.3A, e D = 



0.256!,, 



fk 



T S p(fc) ' 



(A13) 

(A14) 

(A15) 
(A16) 
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where I\(z) is defined in Eq. i|31|) . One can see from Eq. ((A 14(1 that the dominant scattering process is the absorption 
of a phonon with the characteristic energy ~ T as this energy range provides the major contribution to the integral. 

If we neglect the second term of I|A13J| as well as the corrections to the collision integral coming from SCg the 
function f\ can be found right away. For example, if T <C Oz>: 



fl — 



f v 4 



d 



\„4T 4 |fc| dk 



fk(T)\d x T. 



(A17) 



This is the expression for odd (in k) part of the spin boson distribution function. Using it we estimate the omitted 
term of l(A13(l and show that this term is small: 



- A I ? 



Co 



2tt 



kk'(k - k'f 



e (|fe|-|fe'|)/fer _ 1 



-\k'\/ki 



-\k\/k T 



dk' 



Co 



2tt 



kk' 



o(\k\-\k'\)/k T _ 1 



[k 2 - 2kk' + k' 2 ] 



-lfe'1/fex _ i i 
-\k\/k T Z i 



Of three terms in the square brackets only the second gives non-vanishing contribution to the integral. Two others 
do not contribute because they are even in k' and function fl, is odd. Therefore: 



kk'{k-k'f 



|fe|-|/c'|)/fe r _ i 



-\k'\/k T _ i 

-A' 



-\k\/k 2 



2k 
ak^ 



Q / v rlV r) 



e/v 



2k 



T2- lo S 



kn 



d x T,{A18) 



where k min = min{fc<r; Qjj/v}. Since the integral in (|A18(1 diverges at small |fc'| we must introduce small cut-off 
parameter k c . Such a cut-off is associated with the impurity scattering. It is equal to the momentum at which the 
impurity scattering rate becomes comparable with t^, 1 . For now it is enough to assume that k c is not too small so we 

can neglect the logarithm in the above formula. The expression (|A18J) is to be compared with (jA14|) . For T <C &d- 



f 1 1 

i*- a -d x T » 



irak 2 -, 



k'r 



log \ d x T <S4> \k\ 2 < k 2 T . 



(A19) 



In other words, at low temperatures the second term of ((A 13(1 is small as long as \k\ <C kx- The same is true for larger 
temperature (T > 0_d) as well. 

We also need to estimate the integrals over SCg to demonstrate that they are small. First, let us consider the integral 
over SCt. This interval is localized near k! ~ — k: 



SC t = [-(1 + 2c t /v)k, -(1 + 2c t /v)k] U [-(1 - 2c t /v)k, -(1 - 2ci/v)k]. 
It is non- vanishing as long as c; > Ct ■ It is necessary to perform the integral l(A8(l with I = t over SCt '■ 



dk' 

,n.: t 2vr 



kk'(k - k') 



l\2 



^A t {fl-f l _ k ) 
2A t 



e (\k\-\k'\)/k T _ i 

dk' 
2n 



f e W/k T i e -W\,k T _ i 
| el^l/^ - 1 e~\ k \/ kT - 1 



<5C, 



A2 



c 2 t (k~k') 

) 2 (\k\ - \k'\) 2 

U'\\2 _ „2 



1 



e (|fe|-|fe'|)/fc T _ ! 



rfA- 



/U 2 ( | fc| _| fc /| )2 _ c 2 (A ._ A ./ ) 



/\2 



/tfc , W 3 ||fe|- |A;'|| 3 

The absolute value of the last expression can be bounded from above by the following: 



Awe 



4A t 



T\k\ \fl - f_ k \ / dk' (v 2 (\k\ - \k'\) 2 -c 2 (k- k') 2 ) » ^k T k A \fl fl 



1(Q-C t ) 3 (Cl-Ctf 



c 2 



(A20) 



(A21) 



(A22) 



The latter expression can be estimated as oc Atkrk A \fl. |. This quantity is smaller than ((A 14(1 for \k\ < kx as long as 
{k/k T ) 2 <^l. 

Second, we estimate the integral ((A8() with i = I over SCi : 

5Ci = [\k\ + Q Dl /v, \k\ + Q m /v] U [-|*| - Qot/v, -\k\ - Q D i/v\. (A23) 

It is clear without extensive calculations that both terms of l(A8(l are proportional to ($ki) x ~ ( c / v ) 2 ^ 1- Thus, the 
integral over this interval does not contribute significantly to the collision integral of the spin bosons. 

Therefore, we can conclude that for small T the collision integral can be approximated by the expression ((A14(l . 
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APPENDIX B: IMPURITY SCATTERING 



Now we can start evaluating the self-energy corrections to the spin boson Green's function. Our first step is the 
calculation of the Matsubara Green's function T>(t). Since the perturbation is an exponential of the bosonic field <t 
rather than a polynomial we present our calculation in details. The second order correction to the Green's function 
is given by: 

V k (r)-V ok (r)^^ [ dr'dr" ((^(r) + 6_ fe (r)) e-^ 4 ^'-%^*^^ ( 6fe; ( ) + &t_ fe/(0 ))) _ . (Bl) 

nected 

The triangular brackets stand for the Matsubara time-ordered averaging, and xq is the impurity position. The bare 
boson Green's function 



'o*(r)=((&l(T)+6_ fc (r)) (b k (0) + bl k (0))) (B2) 



is a Fourier transform of T>ok,iu — 2wfc/ (cu 2 + ujf.). To calculate the expression in the triangular brackets in Eq. I|B1|) 
we must expand both exponents into Taylor series 



E ( 12 ((to * w ta,^* m (*o.o {km + „_ ( B3 ) 

n,m nected 

and apply Wick's theorem. The external operators b, can be contracted with fields $ in four possible ways 
schematically shown in Fig. (a) both external operators contracted with the monomial $"(xo,t'), (b) both 
external operators contracted with the monomial $ m (:ro, r"), and (c) and (d) one external operator contracts with 
& n (xo, t') while another contracts with $ m (xo, t"). Each contraction with the external operator acts on a monomial 
as a derivative with respect to the field 3> times e %hxa T>Q k j L\k\ or e~ lk x °T>o k i / \/L\k'\. Effectively, contraction with 
an external operator is equivalent to differentiation with respect to $ times some c-number. Because of that one can 
derive: 



o i(k—k')x 



V ok {T')V ok ,(T - T')+V ok (r")V ok ,(T - t") - V ok (T')V ok ,(T - t") - V ok (r")V ok ,(T - r') \ (B4) 



N , e -iV2K<i>(xo,T') e iV2iri(x ,T''] 



The first term in the curly brackets corresponds to the diagram (a) in Fig. [3 the second corresponds to the diagram 
(b) and so on. The gray bubbles in Fig. [2] correspond to ( e -*v^4(o,T') e »V2^*(o,-r'')}. 

As we see from (|B4J| . the impurity scattering explicitly violates the momentum conservation. It will be restored once 
we average over the impurity positions. Such impurity averaging leads to a simple change in the above expression: 
= nS ktk '/a, where n is the dimensionless impurity concentration. After that the impurity position 
xq is arbitrary and can be put to zero. This yields: 

V k (r) - V ok (r) a i J dr' dr" S.V ok (r - T')V Qk {r") + V Qk (r - r")Pofe(r') 

- 2V 0k (T - T")V ok (T")^ k (r' - r"), (B5) 

S fc (r' - r") = — ^ / e -iV2¥i(0,r') e iv^F4(0,r")\ _ (B6) 
air\k\ \ I 

With the help of the standard formula for free bosonic field: 

/ e ai(x,r) e -ai(x',r')\ _ £ a 2 ((*(x,t)-*(x' ,t')) 2 > /2 

the above equation can be rewritten in the w-space as: 

k,iu — A)fc,iw ~ 'E > 0k,iw^'k,wi (B8) 
n 



6J? 



imp 



a7r|A;| 



dr {e™ r - 1) e — ((*(0,r)-*(0,0)) 2 ) _ (B9) 
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The average of the bosonic fields can be evaluated as follows 



g(r) = -( ($(0,t)-$(0,0)) 2 \ = 2 ( $(0, t)$(0, 0) ) - 2 ( $(0, 0) 2 ) (Bid) 



f p iuiT _ 1 



-i /*+oo — a\uj\x/v 

with I(u) = — - / —dx. (Bll) 



7r|u;| Jo 1 + x2 

It is possible to evaluate the integral I(ui) for small (\ui\ <C J ~ u/a) and large (|cj| » J ~ v/a) Matsubara frequency 

m = { 1 'M f or M<<<{ (Bi2) 

v ; [ v/iruj a for \u\ > J. v ; 

The fact that at large \uj\ function I(ui) decays as an inverse u> 2 guarantees the convergence of the Matsubara sum. 
Unfortunately, it is impossible to continue our derivation analytically due to complex properties of /. To overcome 
this difficulty we utilize a trick. Rather than using the actual form of I{oj) we put a different function in its place: 

= -—- . (B13) 

Such a replacement preserves \oj\ <C J properties of the sum but radically improves convergence at large lu. Since we 
are not interested in the high energy (~ J), short time (~ 1/J) properties of £ this procedure is justified. Now the 
Matsubara summation can be done and it leads to: 



iur—\u\/J _ „-\ui\/J 



1 (1 - e - 27rT / J ) 



2 



9(t) T 5 M 2^ l0g (l + e -4-T/./_2e-Wcos2^Tr)' (B14) 

In deriving this formula we used the identity: 

Et-s"!'-'")' (B15) 

UJ>0 

which is correct for any complex 7, Re 7 < 0. To prove this formula one can expand the logarithm in Taylor series 
with respect to powers of e 27rT ' y . 

Once g(r) is calculated it must be substituted into (|B9|) . Now we have to evaluate the Fourier integral: 

F(uj) = / dr - [ ===, ' —■ (B16) 
Jo Vl + e~^ T l J - 2e~ 2 * T ' J cos2ttTt 

This integral can be transformed into an integral over a unit circle in the complex plane: 

F = ^-idz , Z " l ~ 1 (B17) 
2mTJ ^z{{2 + s)z~z 2 -l) 

s = 4tt 2 T 2 /J 2 < 1, z = e 2mTT , m=— integer. (B18) 

2nT 

The polynomial under the square root has three zeros: z = and z 12 ~ 1 ± yfs. Our branch-cut consists of two 
components: the first component stretches from to 1 — y^i, the second component stretches from 1 + ^fs to +00. 
Such a choice of the branch-cut ensures that the integration path, the unit circle, does not cross the branch-cut, see 
Fig. 1101 It is possible now to change the integration path. The new integration path runs counterclockwise along the 
sides of (0,1 — \fs) component of the branch-cut. In turn, this integral can be transformed into an integral over the 
real axis: 



ttT 



F =^ dx ) =. (B19) 



y/(x — l) 2 — SX 
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FIG. 10: The transformation of the integration path in the complex plane for Eq. IB 1711 is shown. The thick gray lines represent 
the branch cuts. The unit circle is transformed into a new contour circumventing one of the branch cuts. This new contour is 
shown by the dash line. 



Finally, the Matsubara self-energy equals to 

" X ^ ,-1 ■ (B20) 

To obtain the retarded self-energy T, R the analytical continuation to — > —ioj/2'KT must be done. For w < T it gives: 

_ z W|p™ (B21) 
^ 2aJ|fc|T V y 

Once the self-energy is calculated one can find the retarded Green's function: 

and extract the lifetime from the position of the Green's function pole: 

u? = u 2 k - 2u; k ^ u . (B23) 

As a result the lifetime T is given by: 

A 2 

r = jr, (B24) 

with the auxiliary parameter A is defined as 

A 2 cx n6J? mp (B25) 
where n is the dimensionless impurity concentration. 
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